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On ' Abstract. We study the limit as A'' ^ oo of the correlations between simultaneous zeros of 

2^ ! random sections of the powers of a positive holomorphic line bundle L over a compact 

■ complex manifold M, when distances are rescaled so that the average density of zeros is 

^ \ independent of N . We show that the limit correlation is independent of the line bundle 

O • and depends only on the dimension of M and the codimension of the zero sets. We also 

, provide some explicit formulas for pair correlations. In particular, we provide an alternate 
derivation of Hannay's limit pair correlation function for SU(2) polynomials, and we show 

, that this correlation function holds for all compact Riemann surfaces. 

S ■ Introduction 

^ \ This paper is concerned with the local statistics of the simultaneous zeros of k random 

O ' holomorphic sections Si, . . . , G H^{M, L^) of the iV*^ power of a positive Hermitian 
0^ . holomorphic line bundle (L, h) over a compact Kahler manifold M (where k < m = dimM). 

The terms 'random' and 'statistics' are with respect to a natural Gaussian probability mea- 
sure du^ on H^{M, L^) which we define below. In the special case where M = CP*" and L is 
the hyperplane section bundle 0{1), sections of correspond to holomorphic polynomials 
of degree A^, and {H^{C¥^, 0{N)), duj^j) is known as the ensemble of SU(m + 1) polynomials 
in the physics literature. To obtain local statistics, we expand a ball U around a given point 
by a factor \/N so that the average density of simultaneous scaled zeros is independent of 
N. We then ask whether the simultaneous scaled zeros behave as if thrown independently 
^ ■ in \fNU or how they are correlated. Correlations between (unsealed) zeros are measured 
^ ! by the so-called n-point zero correlation function K^i^{z^, . . . , z^), and those between scaled 
zeros are measured by the scaled correlation function K'^(-^, . . . , -^). Our main result is 

that the large limits of the scaled n-point correlation functions K^^{-^, . . . , ^) exist 

and are universal, i.e. are independent of M, L and h as well as the point z^. Moreover, 
the scaling limit correlation functions can be calculated explicitly. We find that the limit 
correlations are short range, i.e. that simultaneous scaled zeros behave quite independently 
for large distances. On the other hand, nearby zeros exhibit some degree of repulsion. 

To state our problems and results more precisely, we begin with provisional definitions of 
the correlation functions K^i.{z^, . . . , z^) and of the scaling limit. (See §§0-@ for the complete 
definitions and notation.) In order to provide a standard yardstick for our universality results, 
we give M the Kahler metric u given by the (positive) curvature form of h. The metrics h 
and u then induce a Hilbert space inner product on the space H^{M,L^) of holomorphic 
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sections of L^, for each > 1. In the spirit of |^Z]| we use this £^-norm to define a 
Gaussian probabihty measure du^ on H^{M,L'^). When we speak of a random section, 
we mean a section drawn at random from this ensemble. More generally, we can draw k 
sections (si, . . . , s^) independently and at random from this ensemble. Let ^(si,...,^^,) denote 
their simultaneous zero set and let g^.)! denote the "delta measure" with support on 



^(si,...,Sfc) and with density given by the natural Riemannian volume (2m — 2A;)-form defined 



by the metric u. To define the ra-point zero correlation measure K^i^{z^, . . . , z") we form the 



product measure 

|^(.i,...,.,)r=(|^(si,...,.,)|x---x|Z(,,,...,,,)|) on M":=Mx---xM. 



To avoid trivial self-correlations, we puncture out the generalized diagonal in M" to get the 
punctured product space 

M„ = . . . , 2") G M" : ^ 2« for p q} . 

We then restrict \Z(^si,...,Sk)\"' to M„ and define K^i^{z^, . . . , z"^) to be the expected value 
i?(|Z(sj,...,sj.)|") of this measure with respect to ujy. When k = m, the simultaneous zeros 
almost surely form a discrete set of points and so this case is perhaps the most vivid. 
Roughly speaking, K^f,{z^, . . . , z") gives the probability density of finding simultaneous zeros 



at (^1,...,^") 



The first correlation function Ki^n just gives the expected distribution of simultaneous 
zeros of k sections. In a previous paper by two of the authors, it was shown (among 



other things) that the expected distribution of zeros is asymptotically uniform; i.e. 

K^,{z') = CmkN' + 0{N'-'), 

for any positive line bundle (see ||SZ| , Prop. 4.4]). The question then arises of determining the 
higher correlation functions. As was first observed by |[BBL|| and ||Han|] for SU(2) polynomials 



and by |[BD|| for real polynomials in one variable, the zeros of a random polynomial are non- 



trivially correlated, i.e. the zeros are not thrown down like independent points. We will 
prove the same for all SU(m-|- 1) polynomials and hence, by universality of the scaling limit, 
for any M, L, h. 

To introduce the scaling limit, let us return to the case k = m where the simultaneous 
zeros form a discrete set of points. Since an m-tuple of sections of will have A^"^ times 
as many zeros as m-tuples of sections of L, it is natural to expand f/ by a factor of -\/iV 
to get a density of zeros that is independent of A^. That is, we choose coordinates {zg} 
for which z° = and uj{z^) = ^J^qd-Zq A dzg and then rescale z i— Were the zeros 
thrown independently and at random on U, the conditional probability density of finding a 
simultaneous zero at a point w given a zero at z would be a constant independent of {z,w). 
Non-trivial correlations (for any codimension A; G {1, . . . , m}) are measured by the difference 
between 1 and the (normalized) n-point scaling limit zero correlation function 

K^,^{z\ ...,z-)= hm (c„,Ar^)-" K^J ), (^\ . . . , z") g [/„ . 



Our main result (Theorem p.4| ) is universality of the scaling limit correlation functions: 
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The n-point scaling limit zero correlation function K^^{z^, . . . , z") is given by a 
universal rational function, homogeneous of degree 0, in the values of the function 
^t^{z-w)-^\z-w\ g^^^ j^^g j^^g^ second derivatives at the points {z,w) = (z^^zp'), 

1 < p,p' < n. Alternately it is a rational function in z^, e^'''^'' 

The function e*^(^-'^)-2l^-"'l' which appears in the universal scahng hmit is (up to a constant 
factor) the Szego kernel n^(z, w) of level one for the reduced Heisenberg group H^^^ (cf. §|I|). 
Its appearance here owes to the fact that the correlation functions can be expressed in 
terms of the Szego kernels njv(x, ?/) of . I.e., let X denote the circle bundle over M 
consisting of unit vectors in L*; then Ilisf{x,y) is the kernel of the orthogonal projection 
Un : C\X) n%{X) ^ H%M, L^). Indeed we have (Theorem |2j): 



The n-point correlation K^i^{z^ , . . . , z^) is given by the above universal rational 
function, applied this time to the values of the Szego kernel Un and its first and 
second derivatives at the points {z^, z^ ). 

In view of this relation between the correlation functions and the Szego kernel, it suffices 



for the proof of the universality theorem p.4| to determine the scaling limit of the Szego kernel 



IItv and to show its universality. Indeed we shall show (Theorem p.lD that: 



Let {zi, . . . , Zm, 0) denote local coordinates in a neighborhood U ^ U x of a point 
(z°, X) E X (where {zi, . . . , Zm) are the above local coordinates about zq E M). We 
then have 

The fact that the correlation functions can be expressed in terms of the Szego kernel 
may be explained in (at least) two ways. The first is that the correlation functions may be 
expressed in terms of the joint probability density D^i^{x,^; z)dxd^ of the (vector-valued) 
random variable 

{x,o = [^',e],^,^^, x^ = {s,{zn,...,s,{zn), e = {^s,{zn,...,vs,{zn) 

given by the values of the k sections and of their covariant derivatives at the n points {z^}- 
Our method of computing the correlation functions is based on the following probabilistic 



formula (Theorem 



For N sufficiently large so that the density D^i^{x,^; z) is given by a continuous 
function, we have 



» n 



where ^ = {C,^, . . . and C,^* : L^p — > Tm,zp denotes the adjoint to : Tm,zp 

tN 
^zP- 
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This formula, which is vahd in a more general setting, is based on the approach of Kac 
and Rice |R| (see also [[EK|| ) for zeros of functions on M^, and of ||Hal|| for zeros of 
(real) Gaussian vector fields. Since our probability measure dv^ (on the space of sections) 
is Gaussian, it follows that is also a Gaussian density. It will be proved in § p73| that the 
covariance matrix of this Gaussian may be expressed entirely in terms of IItv and its covariant 
derivatives. This type of formula for the correlation function of zeros was previously used in 
anil and the works cited above. We believe that this formula will have interesting 



applications in geometry. 

A second link between correlation functions and Szego kernels is given by the Poincare- 
Lelong formula. In fact, this was our original approach to computing the correlation functions 
in the codimension 1 case. For the sake of brevity, we will not discuss this approach here; 



instead we refer the reader to our companion article [[BSZ 



From the universality of our answers, it follows that the scaling limit pair correlation 
functions depend only on the distance between points: 



where Kkm depends only on the dimension m of M and the codimension k of the zero set. 
In we give explicit formulas for the limit pair correlation functions Kkm in some special 
cases. Our calculation uses the Heisenberg model, which (although noncompact) is the most 
natural one since the scaled Szego kernels are all equal to Hi , and there is no need in this case 
to take a limit. We also discuss the hyperplane section bundle 0{1) —>■ CP"^, which is the 
most studied, since the sections of its powers are the SU(m + 1) polynomials — homogeneous 
polynomials in m + 1 variables — and the case m = 1 (the SU(2) polynomials) appears 
frequently in the physics literature (e.g., | |BBL| , [b H] , pian| , |K1V1W| , [P'l| ). We give expressions 
for the zero correlations K^j^ for the SU(m + 1) polynomials and by letting ^ cxd, we 
obtain an alternate derivation of our universal formula for the scaling limit correlation. 

We show (Theorem |4.1| ) that Kkm{i^) = 1 + 0(r^e~^ ) as r ^ +oo, and hence these 
correlations are short range in that they differ from the case of independent random points 
by an exponentially decaying term. We observe that when dimM = 1, there is a strong 
repulsion between nearby zeros in the sense that Kii(r) — >■ as r — 0, as was noted by 
Hannay [[Han]| and Bogomolny-Bohigas-Leboeuf ||BBL|| for the case of SU(2) polynomials. 
These asymptotics are illustrated by the following graph (see also ||Han|| ): 




Figure 1. The 1-dimensional limit pair correlation function kh 
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For dimM = 2, the simultaneous scaled zeros of a random pair (si,S2) of sections still 
exhibit a mild repulsion {limr-^o K22{r) = |), as illustrated in Figure 2 below. 
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Figure 2. The limit pair correlation function K22 

The function Kmm{f) can be interpreted as the normalized conditional probability of finding 
a zero near a point given that there is a zero at a second point a scaled distance r from 
(in the case of discrete zeros in m dimensions) . The above graphs show that for dimensions 
1 and 2. there is a unique scaled distance where this probability is maximized. It would 
be interesting to explore the dependence of the correlations on the dimension. To ask one 
concrete question, do the simultaneous scaled zeros in the point case become more and more 
independent in the sense that «:mm(^) — > 1 as the dimension m — > 00? 

When k < m, the zero sets are subvarieties of positive dimension m — k; in this case the 
expected volume of the zero set in a small spherical shell of radius r and thickness e about 
a point in the zero set must be ~ ^^2m-2fc-i^ Hence we have Kkm{i^) ~ r~^'^, for small r. The 
graph of the limit correlation function for the case m = 2, /c = 1 is given in Figure 3 below. 
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Figure 3. The hmit pair correlation function K12 

To end this introduction, we would like to link our methods and results at least heuristi- 
cally to a long tradition of (largely heuristic) results on universality and scaling in statistical 
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mechanics (cf. |FFS|). One may view the rescahng transformation on U as generating 



a renormahzation group. The intuitive picture in statistical mechanics is that the renor- 
mahzation group should carry a given system (read "L M") to the fixed point of the 
renormahzation group, i.e. to the scale invariant situation. We observe that the local rescal- 
ing of U is nothing other than the Heisenberg dilations 6^ on HJ!^^. Since these dilations 
are automorphisms of the (unreduced) Heisenberg group, the Szego kernel of is invariant 
under these dilations; i.e., it is the fixed point of the renormahzation group. As predicted by 
this intuitive picture, we find that in the scaling limit all the invariants of the line bundle, in 
particular its zero-point correlation functions, are drawn to their values for the fixed point 
system (read "Heisenberg model"). 



1. Notation 

We begin with some notation and basic properties of sections of holomorphic line bundles, 
their zero sets, Szego kernels, and Gaussian measures. We also provide two examples that 
will serve as model cases for studying correlations of zeros of sections of line bundles in the 
high power limit. 

1.1. Sections of holomorphic line bundles. In this section, we introduce the basic com- 
plex analytic objects: holomorphic sections and the currents of integration over their zero 
sets. We also introduce Gaussian probability measures on spaces of holomorphic sections. 
For background in complex geometry, we refer to ||GH|| . 

Let M be a compact complex manifold and let L — M be a holomorphic line bundle with 
a smooth Hermitian metric h; its curvature 2-form Qh is given locally by 

(1) e, = -dd\og\\eL\\l, 

where cl denotes a local holomorphic frame (= nonvanishing section) of L over an open set 
U C M, and \\eL\\h = h{eL,eLy/'^ denotes the /i-norm of cl. We say that {L,h) is positive 
if the (real) 2-form to = ^^^^Qh is positive, i.e., if is a Kahler form. We henceforth assume 
that (L, h) is positive, and we give M the Hermitian metric corresponding to the Kahler 
form u} and the induced Riemannian volume form 

(2) dVM = -^u;^ . 

ml 

Since -u is a de Rham representative of the Chern class ci(L) G iJ^(M, M), the volume of 
M equals ^ci(L)'". 

The space H^{M,L^) of global holomorphic sections of = L ® ■ ■ ■ ® L is a finite 
dimensional complex vector space. (Its dimension, given by the Riemann-Roch formula for 
large A^, grows like N"^. By the Kodaira embedding theorem, the global sections of 
give an embedding into a projective space for ^ 0, and hence M is a projective algebraic 
manifold.) The metric h induces Hermitian metrics on given by ||s®^||/iiv = \\s\\^ . 
We give if°(M, L^) the Hermitian inner product 

(3) {suS2)= [ h''{suS2)dVM (si,S2Gi/°(M,L^)), 
and we write |s| = (s, s)^^'^. 
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We now explain our concept of a "random section." We are interested in expected values 
and correlations of zero sets of fc-tuples of holomorphic sections of powers . Since the 
zeros do not depend on constant factors, we could suppose our sections lie in the unit sphere 
in if°(M, L^) with respect to the Hermitian inner product (0), and we pick random sections 
with respect to the spherical measure. Equivalently, we could suppose that s is a random 
element of the projectivization Pi?°(M, L^). Another equivalent approach is to use Gaussian 
measures on the entire space H'^^M, L'^). We shall use the third approach, since Gaussian 
measures seem the best for calculations. Precisely, we give H^{M, L^) the complex Gaussian 
probability measure 

(4) du^is) = —e-^''^'dc, s = Y,CjS^, 



IT' 



where {Sj^} is an orthonormal basis for H^{M, L^) and dc is 2(i7v-dimensional Lebesgue 
measure. This Gaussian is characterized by the property that the 2dN real variables 3ftcj, Qcj 
{j = 1, . . . , d^) are independent random variables with mean and variance |; i.e., 

E Cj = 0, E CjCk = 0, E CjCk = Sjk ■ 

Here and throughout this paper, E denotes expectation. 

In general, a complex Gaussian measure (with mean 0) on a finite dimensional complex 
vector space \^ is a measure u of the form (^, where the Cj are the coordinates with respect to 
some basis. Explicitly, the complex Gaussian measures on C™' are the probability measures 
of the form 

-{A-^z,z) 

(5) ^ dz 

^ ' TT'^detA 

where A = (A-^) is a positive definite Hermitian matrix and 



(C,^) = C-^ = EC 

9=1 



denotes the standard Hermitian inner product in C"^. For the Gaussian measure (^, we 
have 



(6) E {z^Zk) = 0, E {zfzk) = A 



If 1/ is a complex Gaussian on V and t : V ^ V is a surjective linear transformation, then 
t^kZ/ is a complex Gaussian on V. In particular, if = C"^, then, t^^u is of the form (|5|), where 
the covariance matrix A is given by with Zj = Zj o t : V ^ C. 

We shall consider the space S = H'^{M, L'^)^ {1 < k < m) with the probability measure 
dfi = du X ■ ■ ■ X du, which is also Gaussian. Picking a random element of S means picking 
k sections of H^{M, L^) independently and at random. For s = (si, . . . , Sjt) G 5, we let 

Zs = {z e M : si{z) = ■ ■ ■ = Sk{z) = 

denote the zero set of s. Note that if N is sufficiently large so that is base point free, 
then for /i-a.a. s G 5, we have codim Zg = k. (Indeed, the set of s where codim Zs < k is a. 
proper algebraic subvariety of H^{M, L^Y. In fact, by Bertini's theorem, the Zg are smooth 
submanifolds of complex dimension m — k for almost all s, provided is large enough so 
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that the global sections of give a projective embedding of M, but we do not need this 
fact here.) For these s, we let \Zs\ denote Riemannian (2m — 2A;)-volume along the regular 
points of Zg, regarded as a measure on M: 

(7) {\Zsl(p)= I (pdVol2m-2k = — ^—r^ I '^^ 



m- 



rcg 



It was shown by Lelong (see also ||GH|| ) that the integral in (0) converges. (In fact, \Zs\ 
can be regarded as the total variation measure of the closed current of integration over Zg.) 
We regard \Zs\ as a measure- valued random variable on the probability space {S,dfi); i.e., 
for each test function ip e C°(M), (|Zs !,</?) is a complex- valued random variable. 

1.2. Szego kernels. As in |]Z|, we now lift the analysis of holomorphic sections over M 
to a certain bundle X —>■ M. This is a useful approach to the asymptotics of powers of 



line bundles and goes back at least to |BG| 



We let L* denote the dual line bundle to L, and we consider the circle bundle X = {X E 
L* : \\X\\h* = 1}; where h* is the norm on L* dual to h. Let tt : X — > M denote the bundle 
map; if f G L^, then \\v\\h = |(A,f)|, X E X^ = 7r~^{z). Note that X is the boundary of the 
disc bundle D = {X E L* : p(A) > 0}, where p(A) = 1 — ||A||^*. The disc bundle D is strictly 
pseudoconvex in L*, since Qh is positive, and hence X inherits the structure of a strictly 
pseudoconvex CR manifold. Associated to X is the contact form a = —idp\x = idp\x- We 
also give X the volume form 

(8) dVx = -^a A (da)'^ = aA Tx*dVM ■ 

ml 

The setting for our analysis of the Szego kernel is the Hardy space TC'^{X) C C^{X) of 
square-integrable CR functions on X, i.e., functions that are annihilated by the Cauchy- 
Riemann operator Bb (see pp. 592-594]) and are with respect to the inner product 

(9) (Fi, F2) = ^ / F{F2dVx , Fi, F2 E C\X) . 

2vr Jx 

Equivalently, 1-L'^{X) is the space of boundary values of holomorphic functions on D that are 
in £^(X). We let rgx = e*^x (x E X) denote the action on X and denote its infinitesimal 
generator by ^. The 5*^ action on X commutes with db] hence li}{X) = ®'^=q'H%{X) 
where 7i^(X) = {F E 1H?{X) : F{rgx) = e^'^^ F{x)}. A section s of L determines an 
equivariant function s on L* by the rule s(A) = (A, s{z)) (A E L*, z E M). It is clear that if 
r G C then s{z, tX) = rs. We henceforth restrict s to X and then the equivariance property 
takes the form s{rgx) = e^^s{x). Similarly, a section sn of determines an equivariant 
function sn on X: put 

sn{X) = {X^^,sn{z)) , AgX,, 

where A*^^ = A® - • - C^A; then SNirex) = e^^^sxix). The map s ^— s is a unitary equivalence 
between if°(M, L^) and nlf{X). (This follows from (|)-(|) and the fact that a = d9 along 
the fibers of vr : X M.) 

We let IItv : 'C^(X) T-C%{X) denote the orthogonal projection. The Szego kernel 
IlN{x,y) is defined by 

(10) njvF(x) = / n^(x, y)F{y)dVx{y) , F E C^{X) . 

Jx 
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It can be given as 



(11) u^{x,y) = J2s^i^)sf{y): 

where , . . . , S^^ form an orthonormal basis of H^{M^L^). Pick a local liolomorphic 
frame for L over an open subset U C M, let denote the dual frame, and write 
h{z) = h{eL{z),eL{z)) = IIclH^. The map {z,e^^) i— >• e^^h{zy^'^e*];^{z) gives an isomorphism 
U X ^ 7r~^{U) C X, and we use the coordinates {z, 9) to identify points of 7r~^(C/). For 
s e H^{M, L^), we have 

(12) s{z, d) = {s{z), e'^^'hizf/^eliz)) = e'^^hiz)""^ f{z), s = /ef ^ . 

Although the Szego kernel is defined on X, its absolute value is well-defined on M as 
follows: writing — f^ef^, we have 

(13) Ur,{z,e;w,^) = e^^('^-^)n^(^,0;^/;,0) = e^''^'-'^^h{zf/'h{wf/'Y.f^^^^'WM^ 

for z,w ^ U. (Here we may take U to be the disjoint union of connected neighborhoods of 
z and w, if 2; is not close to w.) Thus we can write 

\IIn{z,w)\ = \Un{z,0;w,0)\ , 

which is independent of the choice of local frame cl- On the diagonal we have 

n^(^, z) = Ur,{z, 9; z,9)=J2 ll-^f (^) lU- • 

The Hermitian connection V on L induces the decomposition Tx = (BTx into horizontal 
and vertical components, and we let denote the horizontal lift (to X) of a vector field t 
in M. We consider the horizontal operators on X: 

"Zq ~ "^{d/dZq)" , % - a{d/dZq)H , 

where zi, . . . , Zrm9 denote local coordinates on X. We note that 

(14) <5=(V^.)^, seH\M,L''), 
where is the induced connection on . We then have 
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^ / \ 



(15) dldlUr,{z,e;w,^) = ^ (^S^f )» )^(^) 



e^^(^"-)M.)^/^M^)^/^£/j;(.)/jv 



i=i 

We can also use (1121) and (O) to describe the horizontal lift in local coordinates: 



^ ' dz, 2 dz, 09- 

1.3. Model examples. In two special cases we can work out the Szego kernels and their 
derivatives explicitly, namely for the hyperplane section bundle over CP™ and for the Heisen- 
berg bundle over C™, i.e. the trivial line bundle with curvature equal to the standard sym- 
plectic form on C". These cases will be important after we have proven universality, since 
scaling limits of correlation functions for all line bundles coincide with those of the model 
cases. 

In fact, the two models are locally equivalent in the CR sense. In the case of CP*", the circle 
bundle X is the 2m + 1 sphere 5"^™"*"^, which is the boundary of the unit ball _B2m+2 ^ 
In the case of C"^, the circle bundle is the reduced Heisenberg group H™^, which is a discrete 
quotient of the simply connected Heisenberg group C™ x M. As is well-known, the latter is 
equivalent (in the CR and contact sense) to the boundary of 52m+2 (|g|]). 

1.3.1. SU(m + 1) -polynomials. For our first example, we let M = CP™ and take L to be 
the hyperplane section bundle 0{1). Sections s G iJ°(CP™, (9(1)) are linear functions on 
(^m+i. 2ero divisors Zg are projective hyperplanes. The line bundle 0{1) carries a natural 
metric /ips given by 

(17) PIUps(M) = , W = {wo, . . . , Wrn) G C^+' , 

\w\ 

for s e C™+^* = if°(CP™,0(l)), where |wp = Er=o kjf and [w] e CP™ is the complex 
line through w. The Kahler form on CP™ is the Fubini-Study form 

(18) ujps = ^e,,3 = log \w\' . 

The dual bundle L* = 0{—l) is the affine space C™^^ with the origin blown up, and 
X = 5^™+! c C™+^ The A^-th tensor power of 0(1) is denoted 0{N). Elements sn e 
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H'^[C¥"^,0{N)) are homogeneous polynomials on C™'"^^ of degree A^, and sn = SAr|5'2m-i. 
The monomials 



(19) 



{N + m)\ 



TT^Jo! ■■■Jr. 



= (jO,--- ,Jm), \J\=^ 



form an orthonormal basis for ff°(CP"^, C(A^)). (See [||, §4.2]; the extra factor (^)^^^ 
in ( p!9D comes from the fact that here CP*" has the usual volume whereas in ||SZ|| , the 
volume of CP™ is normalized to be 1.) Hence the Szego kernel for 0{N) is given by 



(20) 

Note that 



^ / N ^ (N + m)\ J J (N + m)\, 



N 



TT'-Jo! ■■■Jr. 



Il{x,y) = IlN{x,y) = — {x,y)) (™-+^) = 27r x [classical Szego kernel on S 

N=l ^ 



(The factor 2tt is due to our normalization 



1.3.2. The Heisenberg model. Our second example is the linear model x C ^ C™ for 
positive line bundles L ^ M over Kahler manifolds and their associated Szego kernels. It is 
most illuminating to consider the associated principal bundle C™ x 5*^ — > C", which may 
be identified with the boundary of the disc bundle D C L* in the dual line bundle. This 
bundle is the reduced Heisenberg group HJ^^ (cf. p. 23). 

Let us recall its definition and properties. We start with the usual (simply connected) 
Heisenberg group H"* (cf. |^| note that different authors differ by factors of 2 and tt 
in various definitions). It is the group x M with group law 

{C,t)-{r],s) = {C + r],t + s + Q{C-r])). 

The identity element is (0,0) and {(,t)~^ = {—(,—t). Abstractly, the Lie algebra of 
is spanned by elements Zi, . . . , Z^, Zi, . . . , Z^, T satisfying the canonical commutation re- 
lations [Zj, Zk] = —iSjkT (all other brackets zero). Below we will select such a basis of left 
invariant vector fields. 

H™ is a strictly convex CR manifold which may be embedded in C'""'"^ as the boundary of 
a strictly pseudoconvex domain, namely the upper half space U"^ := {z E C""*"^ : Qzm+i > 
i Er=i kjf }• The boundary of equals dU"" = {z e C"*+^ : "^Zm+i = \ EJli kjf }• H"^ 
acts simply transitively on dW^ (cf. ||S|], XII), and we get an identification of H*" with dW^ 
by: 

[C,t]-(C,t + ^ICI')eaw™. 

The Szego projector of H*" is the operator H : /^^(H™) 7i^(H'") of orthogonal projection 
onto boundary values of holomorphic functions on which lie in C? . The kernel of H is 
given by (cf. XII §2 (29)) 

Jl{x,y) = Kiy-'x), K{x) = -C^^[t + t\C\Y"^ G P'(H™) . 
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The linear model for the principal bundle described in §1.2 is the so-called reduced 
Heisenberg group H^^ = H"^/{(0, 2Trk) : k e Z} ^ C"" x with group law 

It is the principal 5"^ bundle over C"* associated to the line bundle Lh = x C. The metric 
on Lh with curvature © = dd\z\'^ is given by setting hu{z) — e~l^l ; i.e., \ f\hH — 1/1^"'^' 
The reduced Heisenberg group H]^j may be viewed as the boundary of the dual disc bundle 
D C and hence is a strictly pseudoconvex CR manifold. 

It seems most natural to approach the analysis of the Szego kernels on H™^ from the 
representation-theoretic point of view. Let us begin with the case = 1. We thus consider 
the space Vi C >C^(HJ^^) of functions / satisfying = /, which forms a (reducible) 

representation of H™^ with central character e*^. By the Stone- von Neumann theorem there 
exists a unique (up to equivalence) representation {Vi, pi) with this character and by the 
Plancherel theorem, Vi = Vi ® V{. 

The space of CR functions in Vi is an irreducible invariant subspace. Here, by CR functions 
we mean the functions satisfying the left-invariant Cauchy-Riemann equations Z^f = on 
H™jj. Here, {Z^} denotes a basis of the left-invariant anti-holomorphic vector fields on 
H^^. Let us recall their definition: we first equip H™^ with its left-invariant contact form 
= '^q{uqdvq — VqclUq) + d9 {( = u + iv) . The left- invariant CR holomorphic (resp. 
anti-holomorphic) vector fields Zq (resp. Zq) are the horizontal lifts of the vector fields ^ 
(resp. ^) with respect to a^. They span the left-invariant CR structure of H^^ and the 
Zq obviously have the form Zq = ^ + where the coefficient A is determined by the 
condition a^{Z^) — 0. An easy calculation gives: 

yL _ ^ i - d yL _ ^ ^ ^ 

The vector fields {^, Z^, Z^} span the Lie algebra of and satisfy the canonical com- 
mutation relations above. 

We then define the Hardy space 7i^(H™jj) of CR holomorphic functions, i.e. solutions of 
Z^f = 0, which he in C^iHl'^^)- We also put Hj = Vin7^2(H™d)- The group H;:^^ acts by left 
translation on Tif. The generators of this representation are the right-invariant vector fields 
Zq,Zq together with ^. They are horizontal with respect to the right-invariant contact 
form = 'Y^q{uqdvq — Vqduq) — d6 and are given by: 

yR "9 i _ d d id 

''^diq~ 2^' 06' ^dFq^ 2^' 86 ■ 

In physics terminology, Z^ is known as an annihilation operator and Z^ is a creation oper- 
ator. 

The representation is irreducible and may be identified with the Bargmann-Fock 
space of entire holomorphic functions on C'" which are square integrable relative to e^'^' 
(or equivalently, holomorphic sections of the trivial line bundle Lh = C™ x C mentioned 
above, with hermitian metric h-n = e^'^'^). The identification goes as follows: the function 
(Pq{z, 6) := e*^e~'^' is CR holomorphic and is also the ground state for the right invariant 
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"annihilation operator;" i.e., it satisfies 

Any element F{z,6) of Tif may be written in the form F{z,6) = f{z)(fQ. Then Z^F = 
(-^f)LpQ, so that F is CR if and only if / is holomorphic. Moreover, F E £^(HJ^j) if and 

only if / is square integrable relative to e~'^'^. 

The Szego kernel U^{z,6,w,ip) of H™^ is by definition the orthogonal projection from 
/:(H^d) to Hf. As will be seen below, Uf{z,e,w,<^) = -i-e^(^-'^)e(^-'^-^l^l'-5l"'l'), which is 
the left translate of (po by {—w,—ip). In the physics terminology it is the coherent state 
associated to the phase space point w. 

So far we have set = 1, but the story is very similar for any A^. We define T-C% as the 
space of square- integrable CR functions transforming by e*^^ under the central S^. By the 
Stone-von Neumann theorem there is a unique irreducible Vat with this central character. 
The main difference to the case A^ = 1 is that Tij^ is of multiplicity A^"^. The Szego kernel 
n^(x, y) is the orthogonal projection to T-C% and is given by the dilate of 11^. Thus, 

To prove these formulae for the Szego kernels, we observe that the reduced Szego kernels 
are obtained by projecting the Szego kernel on to as an automorphic kernel, i.e. 

nH(x,y) = ^n(x,y(0,27rn)). 

Let us write x = {z,6),y = {w,!f). Then the A^-th Fourier component U^{x,y) of 11^, i.e. 
the projection onto square integrable holomorphic sections of , is given by: 

U^{x,y) = I e-'^'Tl{(^'x,y)dt = I e-'^'K{e''y-^x)dt 
Jr Jr 



Here we abbreviated the element (0, e**) by e**. Change variables t\-^t — 9 + ip — Q{z ■ w) to 
get 

nH(x, y) = e'N{e-^)^iN^iz-n^^ e-'^'K{z - w, t)dt 

where Kt is the Fourier transform of K with respect to the t variable. By p. 585], the full 
M X R Fourier transform of K is given by k{z, N) = C'^e-\'\'/^^, so by taking the mverse 
Fourier transform in the z variable we get the Fourier transform just in the t variable: 

(21) nH(a;,y) = ^Ar-e^^(^-^)e*^^(^-'^)e-5^l^-"'l'. 

(Our constant factor in (21) is determined by the condition that 11^ is an orthogonal 
projection.) 
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In our study of the correlation functions, we will need explicit formulae for the horizontal 
derivatives of the Szego kernel. The left-invariant derivatives are given by 

(22) 

N-^Z^W^U^iz, e- w, = N\z,. - w,,){w, - z^)U^{z, 9; w, + N5,,,U^{z, 9; w, ^) . 

Comparing the definitions of the horizontal vector fields with (0), using h-^ = e"'^'^, we 
see that = Z^, as expected, since agrees with the contact form a for Lh (as defined 
in §1.2| ). We will see later that our formulas for computing correlations are valid with any 
connection, and thus it is sometimes useful to also consider the right invariant derivatives: 

N-^Z^U^{z,9;w,^) = Nw,U^{z,9;w,^), 

(23) 

N-"'Zj^W^'^U^{z,9;w,^) = N^z,,WgU^{z,9-w,<f) + N6gg>U^{z,9;w,<f). 

Remark: Recall that the metric on 0{N) —>■ CP"* is given by {z) = (1 + |2;p)~^ using 
the coordinates and local frame from Example |1.3.1 . Since 

h''{z/VN)^hn{z), 

the Heisenberg bundle can be regarded as the scaling limit of 0{N). (Of course, in the same 
way Lh is the scaling limit of , for any positive line bundle L — ^ M.) 



2. Correlation functions 

This section begins with a generalization to arbitrary dimension and codimension a formula 
of |[Han|| and ||BD|| for the "correlation density function" in the one-dimensional case. In fact. 



our formula (Theorem ^TTP applies to a general class of probability spaces of /c-tuples of (real 
or complex) functions. We then specialize to the case where the space of sections has a 
Gaussian measure. Finally, we show how the correlations of the zeros of fc-tuples of sections 
of the A^-th power of a holomorphic line bundle are given by a rational function in the Szego 
kernel n^v and its derivatives (Theorem p.4| ). 



2.1. General formula for zero correlations. For our general setting, we let {V,h) be 
a Hermitian holomorphic vector bundle on an m-dimensional Hermitian complex manifold 
{M,g). (Here, we make no curvature assumptions.) Suppose that 5 is a finite dimensional 
subspace of the space H^{M, V) of global holomorphic sections of V, and let dfi be a prob- 
ability measure on S given by a semi-positive C° volume form that is strictly positive in a 
neighborhood of G 5. (We shall later apply our results to the case where V = L^©- ■ -©L^, 
for a holomorphic line bundle L over a compact complex manifold M, and S = H^{M, V) 
with a Gaussian measure dfi. Our formulation involving general vector bundles allows us 
to reduce the study of n-point correlations to the case n = 1, i.e., to expected densities of 
zeros.) 

As in the introduction, we introduce the punctured product 

Mn = {{z\...,z'') e M X ■■■ X M : zP ^ z'^ forp^q}, 
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and we write 

s{z) = {s{z'), s{z^)) , Vs{z) = {Vs{z'), Ws{z-)) , z = {z\ . . . , z") e M„ , 

where Vs(C) £ (8) is the covariant derivative with respect to the Hermitian connection 
on V. We define the map 

J:M„ [(CerXf)®'^]", J{z,s)^{s{z),Vs{z)); 

i.e., J^{z, s) is the 1-jet of s at 2; G Af„. 

We write g = Qqq'^Zq ® dzq>, hjj> = h{ej, Cj), where {zi, . . . , z^} are local coordinates 
in M and {ei, . . . , e^} is a local frame in (m = dim M, k = rank V). We let G = det{gqgi), 
H = det(/ijj/). We let 

dC = —io^ = G(C) n d^QdQQ , CeM 

i=i 

denote Riemannian volume in M, and we write 

(24) x^^Yl ^j^i ^^^'> ' ^ ^^^^'> n dMf^d'^^ e , 

j j 

= E "i'?^^^ ® ^^-1-^' '^^^ = G(zf )-'=i/(z^')- JI dKa^^d^a^^ e,- e (71/ ® ^).^ ■ 

The quantities dx^, d^^ are the intrinsic volume measures on V^p and {TIj^V)zp, respectively, 
induced by the metrics g,h. 

Definition: Suppose that is surjective. We define the n-point density Dn{x,^, z)dxd^dz 
of II by 

(25) J*{dz X d/i) = Dn{x, ^, z)dxd^dz , x = (x\ . . . , x") e x • • • x , 

e = e (t;,®1^).i x---x(T^®y),n, z = (^\...,z") eM„, 

dx^dx^--- dx" , di^di^--- dC , dz^dz^--- dz"" . 
In this case, for each 2; G M„, the (vector- valued) random variable {s{z),'Vs{z)) has (joint) 
probability distribution Dn{x,^, z)dxd^. 

Remark: If we let n = 1 and fix a point z e M, then the measure D{x,^, z)dxd^ is 
intrinsically defined as a measure on the space J}:{M,V) of 1-jets of sections of V at z. 
Taking a section to its 1-jet at z defines a map Jz '■ S ^ Jl{M,V) and hence induces 
a measure JT^*/^. on J^(M, V^) independently of any choices of connections, coordinates or 
metrics. Similarly for n > 1, £)(x, ^, z)dxd^ is an intrinsic measure on 11^=1 Jlp{^-i 

For a vector-valued 1-form ^ e ® = Hom(TM,z, V^), we let ^* e Hom(V2,TM,2) 
denote the adjoint to ^ (i.e., = {v,it) ), and we consider the endomorphism e 

Hom(T4, Vz). In terms of local frames, if 

3 3,1 

then ^ 
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where ['jqq') = {gqq') ^; hence we have 

(26) = hjijiiajq-^qiqajiiqi Cj ® e*, . 



J,j'J",ci,Q' 

Its determinant is given by 

(27) det(er) =HdetlYl ^n^g'.ajW ] = H det{Q, 0') = ^IKi A ■ ■ ■ A e 

^9''?' ^ l<j,j'<k 



2 

fcll 



Remark: The measure det{^^*)D{0,^, z)dC,dz will play a fundamental role in our study of 
correlation functions. We observe here that it depends only on the metric uj on M, and in 
the case where the zero sets are points {k = m), it is independent of the choice of metric on 
M as well. Indeed, as mentioned in the previous remark, D{x,C,, z)dxd^ is well-defined on 
J^(M, V"). The conditional density D{0,^, z)d^ equals Jz*!^'/ dx\x=o and thus depends only 
on the choice of volume forms dx'^ on VzP. Since dz/dx transforms in the opposite way to 
det it follows that det(^^*)Z^(0, ^, z)d^dz is an invariantly defined measure on (T^ ® K)". 

Recall that for s G 5 so that codimZ^ = /c, we let \Zs\ denote Riemannian (2m — 2k)- 
volume along the regular points of Z^, regarded as a measure on M. 

Definition: For s G iS so that codimZg = /c, we consider the product measure on M„, 

\Zs\ = ( \Zs\ X ■ ■ ■ X \Zs\ ) . 

^ V ' 

n 

Its expectation E is called the n-point zero correlation measure. 

We shall use the following general formula to compute the correlations of zeros and to 
show universality of the scaling limit: 

Theorem 2.1. Let M, V, S, dfi be as above, and suppose that J is surjective and the volumes 
\Zs\ are locally uniformly bounded above. Then 

/n 
d^D^iO,^;z)lldet ieCl ■ 
p=i 

The function Kn{z^, . . . , 2;"), which is continuous on M„ is called the n-point zero correla- 
tion function. For k < m, (p8|) holds on all of the n-fold product M x ■ ■ ■ x M, including the 
diagonal, and Kn is locally integrable on M x ■ ■ ■ x M (and is infinite on the diagonal). In 
the case k = m, when the zero sets are discrete, the zero correlation measure on M x ■ ■ ■ x M 
is the sum of the absolutely continuous measure Kn{z)dz plus a measure supported on the 
diagonal. 



Proof of Theorem |^. ]\ : Consider the Hermitian vector bundle Vn = ©p=i T^pV — > Mn, 



where Hp : M„ — > M denotes the projection onto the p-th factor. By replacing V ^ M with 
Vn Mn and s G H^{M, V) with 



..,s(z"))G//°(M„,V; 



siz ,...,z ] = [siz 
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and noting that Tm„,z = Y[p'^M,zP and = \Zs\, we can assume without loss of generahty 
that n = 1. 

It follows from the above remarks that D{0,^;z) does not depend on the choice of con- 
nection on V. We can also verify this in terms of local coordinates: write s = Yl^j^j^ 
Vs = ^ ajqdzq (S> Cj as in (pSD; we have ajg = I7- + J2k ^^dj^- Then if we write cl% = we 
have 

Hence D{0,^; z) is unchanged if we substitute the (local) flat connection given by a^^. 

We now restrict to a coordinate neighborhood U <Z M where V has a local frame {cj}. By 
hypothesis, we can suppose that the Cj are restrictions of sections in S. We write s = ^ SjCj, 
and by the above we may assume that Vs = ^ dsj ® e^. We use the notation 

= Vdet(a*) , for e G T^,, ®V, = Hom(TM,., V,) . 



Then by (|2 



llVsll 



H\\dsi A ■ ■ • A dsk\ 



1^1 



where ^ is the (/c, A; ) -form on U given by: 



^ = if 

Thus, by the Leray formula. 



(29) 



\dsi A ■ ■ ■ A dsk\ 



dsi A dsi ) A ■ • • A ( -dsk A ds 



dz 



^dsi A (isi ■ ■ ■ A ^dsk A c/s^ 
Define the measure A on M x 5 by 



„^ ,,,9 dz 
Vsf — 



(30) 
Then 



7r*A = E IZ'J" 



where tt : M x S M is the projection. Hence, 



(31) 



A 



dfi{s) \Zs 



dfi{s) 



iVsi 



,dz 

¥ 



For (almost all) x G C^, let I{s,x) be the measure on f/ given by 



{I{s,x),ip) 



V9(^)(iVol(2m-2fc)n(2) 



,dz 



where the second equality is by ( P^ applied to s — Yli^j^ji.z). Then 



(32) 



I{s,x)dx = l\Vs{z)fdz. 



Now let Xx be the measure on U given by 



(Ax, V^) = / {I{s,x),(p)dfj.{s) 
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Claim: The map x (A^, ip) is continuous. 

To prove this claim, we first note that the hypothesis that \Zs\ is locally uniformly bounded 
implies that {I{s,x),ip) < C < +00 uniformly in s,x. Thus we can assume without loss of 
generality that fi has compact support in S. By hypothesis, the map 

or : t/ X 5 ^ C^ a{z, s) = {si{z),..., Sk{z)) 

is a submersion. We can now write A^^ as a fiber integral of a compactly supported form: 

(m - A;)! J^-^^^ 

and thus A^ is continuous, verifying the claim. 

We note that Aq = A|;7. Hence, to complete the proof, we must show that 

7r*Ao = Ki{z)dz\u ■ 
By ( p5|) and (|32|) , for a test function (^^(x, ^,2;), 

^{x,i,z)UfDr{x,i,z)dxdidz = [dfi{s) [^{J{z,smVs{z)fdz 



dx J {I{s,x),(p o J')dij{s) 
{Xx,<fo J)dx. 

By choosing {p{x^^,z) = pe{x)ip{z), where ps is an approximate identity, and letting e — > 0, 
we conclude that 



ij{z)K,{z)dz = / ij{z)\\ifD,{Q,i,z)didz = {Xo.^Piz)) 



□ 



We note the following analogous formula for real manifolds: 



Theorem 2.2. Let V he a C°° real vector bundle over a C°° Riemannian manifold M, and 
let p be a probability measure on a finite dimensional vector space S of C°° sections of V 
given by a semi-positive volume form that is strictly positive at 0. Suppose that the volumes 
\Zs\ are locally uniformly bounded above. Let Dn{x,^, z)dxd^dz denote the n-point density 
of p. Then 



(33) E 



\Zs\' = Kn{z)dz, Kn{z)= / deDfc(0,e,^)n Vdet(ee*) 



The proof is similar to that of Theorem 2J, except that (^) is replaced by the Leray 
formula 

(34) |Z,| = ||dsi A ■ ■ ■ A rfsfc" 



dsi A ■ ■ ■ A dsk 



in the real case. 



UNIVERSALITY AND SCALING OF CORRELATIONS BETWEEN ZEROS 19 



2.2. Formula for Gaussian densities. We now specialize our formula from Theorem 2.1 



to the case where /i is a Gaussian measure. Fix z = {z^,. . . , z") G M„ and choose local 
coordinates {z^} and local frames {e^} near z^, p = l,...,n. We consider the random 
variables t^, a^^ given by 

k km 

(35) s{zn = Y.h^,e^,, Vs(z^) = J]J]a^,ciz^®e^, p=l,...,n. 

j = l j = l q=l 

By (§-(1) and (|2|)-(^ the n-point density 



Dn{x, ^, z)dxd^dz = Dn 



WG{z^)-^H{zP) 
Lp=i 



dhdadz 



fcn(l+m)^g^A ' \a 



is given by: 

(36) D^{b,a;z) 
where 

(37) ~ \ Tj* r' 

= (A- ) = {Eh^,), B^= [Bf:,^ = ( E ) , CV. = (C-,) = ( E a^^,, ) ; 

= = l,...,n; g, g' = 1, . . . , m. 

(We note that An-, Bn, Cn are kn x /cn, /cn x fcnm, knm x /cnm matrices, respectively; j, p, g 
index the rows, and g' index the columns.) 

The function -D„(0, a; z) is a Gaussian function, but it is not normalized as a probability 
density. It can be represented as 

(38) D^{Q,a-z) = ZrXz)D^Sa-z), 
where 

(39) D^^X(^-z) = -^--^--^exp(-(A;ia,a)) 
is the Gaussian density with covariance matrix 

(40) K = C^-BlA-^Bn=ic^^X'' E r^^P^^iw 

V ii,pi,i2,P2 / 

and 

det A„ 1 



(r = A-^) 



^^^^ ^n(2) - ^fcn^g^^^ 7r'^«detA„ 



This reduces formula ( PBD to 

1 / " 

\p=l / A 

where (■)a„ stands for averaging with respect to the Gaussian density D\^{a; z), and (7^^/) 
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2.3. Densities and the Szego kernel. We return to our positive Hermitian line bundle 
(L, h) on a compact complex manifold M with Kahler form u = |6/i- We now apply formulas 
to the vector bundle 



k 

and space of sections 

S = H%M, V) = H'^{M,L^)^ 

with the Gaussian measure /i = z/^v x • ■ ■ x z/jv, where z/^v is the standard Gaussian measure 
on H^{M, L'^) given by We denote the resulting n-point density by -D^, and we also 
write A„ = A^^, A„ = A^^., etc. 

As above, we fix z = {z^, . . . , z"') G M„ and choose local coordinates {z^} near z^\ p = 
1, . . . ,n. We also choose local frames {e^} for L near the points z^ so that 

\\eliznh = l. 

For s G 5, we write 

(43) s{z^) =1 : = ; I (ei(z^))«^ 



(44) V^.,(^^) = a^^dzl ® {eliznr'^ . 

9=1 

Since the sj are independent and have identical distributions, we have 

(45) A^, = {Af^^ = (5,,,E (fcrfef)) , B^, = {Bf^,^,) = (5,,,E (6?<,)) 

c^:.=(c:;';v) = fe'E(a?x;'))- 

We write 



a=l \a=l / 



where {5*^} is an orthonormal basis for H^{M,L^). Using the local coordinates {z^,0) in 
X as described in §L1, we have by (|45|) and ([131) (noting that /i(2;^) = by the above choice 
of local frames) 



5 

(In 



(46) A^^, = E {caC(s)mznft{zP') = /«(^')/" (^^') = ^n'^N{z^, 0; z^', 0) . 

a,/3=l a=l 

Similarly, 



(47) = Yl f'i^nCizP') = 5,yd^n^iz^, 0; z^', 0) 



(48) Cjj^^, = Y fU^nitA'^') = 5,y<e,n^(z^ 0; z^' , 0) . 
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Lemma 2.3. There is a positive integer Nq = No{M,n) such that 

det (n^(^^0;^^',0))^^^^^,^^^0, 
for distinct points Z"*^, . . . , of M and for all N > Nq. 

Proof. It is a well-known consequence of the Kodaira Vanishing Theorem (see for example, 
||GH]| ) that we can find A'^o such that if iV > A'o and xi, . . . ,Xn € M with Xp ^ x\ for 
2 < p < n, then there is a section s G H^{M.,L^) with s{xi) 7^ and s{xp) = for 
2 < p < n. 

We write Appi = IIn{zP,0; z^' ,0). Suppose on the contrary that det{Appi) = 0, and chose 
a nonzero vector v = (t>i, . . . , t>„) such that Ylp'^-'p^pp' — 0- Then recalling (p!T|), we have 

\x P 



(49) = = E ^P^" (^'' 0)^« (^"'' 0) V = E 

p,p' P,p',ct a=l 

where Xa = ^p'VpS^{zP^O). Since the span H'^{M, L^), it follows that for all s E 
H^{M, L^), we have 'Yl,p''^p^i.'^^) = 0- But this contradicts the fact that, choosing po with 
fpo 7^ 0, we can find a section s G H^{M, L^) with 5(2;^°) 7^ and s{zP) = for p 7^ po- CH 

Thus we see that the n-point correlation functions depend only on the Szego kernel, as 
follows: 

Theorem 2.4. Let {L, h) he a positive Hermitian line bundle on an m- dimensional compact 
complex manifold M with Kdhler form uj = |0/i, let S = H^{M, L^)^ (k > 1), and give S 
the standard Gaussian measure /i described above. Let n > 1 and suppose that N is suffi- 
ciently large so that J is surjective. Let z = {z^, . . . , z") G Mn and choose local coordinates 
(Ci) • • • ) Cm) (it coc/i point such that Qhiz^) = d(q A dQ, 1 < p < n. Then the n-point 
correlation K^f,{z) is given by a universal rational function, homogeneous of degree 0, in the 
values ofT^N o-nd its first and second derivatives at the points {z^,zP ). Specifically, 

^ _ p.,^(n^(z^z^'),C,n^(^^^^'),<n^(^^^"'),« n^(^^^"')) 

(50) ^„fc(^) - - k{n+l) 

det (n^(^^^^')), 



(1 < p,p' < n, 1 < q, q' <m), where Vnkm is a universal homogeneous polynomial of degree 
kn{n + 1) with integer coefficients depending only on n, k, m. 

Proof. The n-point zero correlation K^j^{z) is given by equation (^2]) with 7^^, = 5qqi. By the 
Wick formula (|^, (1.13)]), the expectation 

' n 

JJdet (aP*aP) 

in ( ^2|) is a homogeneous polynomial (over Z) of degree kn in the coefficients of A„. By ( ^0|) 
and (^), the coefficients of det ijlj^{z^ , z^ ))A„ are homogeneous polynomials of degree n + 1 
in the coefficients of An, Bn, Cn- The conclusion then follows from (146|)-(^8D. □ 
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Remark: In the statement of Theorem 2A, we wrote IIn{z, w) for IIn{z, 6; w, (f). Since the 
expression is homogeneous of degree 0, it is independent of 9 and (p. Alternately, we could 
regard IIn{z, w) as functions on M x M having values in ® (replacing the horizontal 
derivatives with the corresponding covariant derivatives); again the degree homogeneity 
makes the expression a scalar. Furthermore, since Theorem ^TT] is valid for all connections, 
we can replace the horizontal derivatives in (^) with the derivatives with respect to an 
arbitrary connection. 

2.4. Zero correlation for SU(m + l)-polynomials. In this section, we use our methods 
to describe the zero correlation functions for SU(m + l)-polynomials. We do not carry out 
the computations in complete detail, since we are primarily interested in the scaling limits, 
which we shall compute in 

The SU(m + 1) -polynomials are random homogeneous polynomials of degree iV > on 

(51) s{z) = s{zo,zu...,Zm) = Yl VN\/J\cjz\ z^ = z^' ■ ■ ■ z^ , J! = Jo!---Jm!, 

|J|=JV 

where the coefficients cj are complex independent Gaussian random variables with mean 
and variance 1: 

(52) Ecj = 0; Ecjc]^ = 5jK, ^jk = ^jako ■ ■ ■^jn.km] EcjC/^ = 0. 

Then s{z) is a Gaussian random polynomial on C^^^ with first and second moments given 
by 

(53) Es(z) = 0; E s{z)s{w) = {z, w)^ = I ^ ZgW^ j ; Es{z)s{w) = 0. 

This implies that the probability distribution of s{z) is invariant with respect to the map 
s{z) s{Uz) for all U e SU(m + 1). 

Let (iSAr,/iAr) denote the Gaussian probability space of independent /c-tuples {k < m) of 
SU(m + l)-polynomials of degree N . For s = (si, . . . , s^) G Sjy, the zero set 

Zs = {z : si{z) = ■■■ = Skiz) = 0} . 

is an algebraic variety in the complex projective space CP™. We will assume that CP™" is 
supplied with the Fubini-Study Hermitian metric u, which is SU(m + l)-invariant. In the 
affine coordinates z = {1, Zi, . . . , z^), 

J2dZqAdz^ (E '^dzg) A {Y, ZgdZq) 



(54) uj = \^dd log (l + Yl 



|2^2 



I.e., 



(55) UJ= ^^^^ ^ gqq'dZg A dZq, , 



(1 + 


z 




' - ZqZqi 


(1 + 


z 


2^2 



To simplify our computations, we consider only points z^ with finite affine coordinates, 
z^ = (1, Zi, . . . , z^), p = 1, . . . ,n, and we regard the SU(m+l)-polynomials sj as polynomials 
of degree < on C™; i.e., we regard the sj as sections of the trivial line bundle on C" with 
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the flat metric h = 1 (so that the covariant derivatives coincide with the usual derivatives of 
functions). 

As above, we consider the random variables 



and we denote their joint distribution by 

D^,{b,a;z)dbda, 6 = . . . , 6") , 6^ = (6p,=i,...,fc; 

(56) 

a = (ai,...,a"), = • 

(Here, the n-point density is with respect to Lebesgue measure db = Yl dt^dhj, da = 
Ylda^gdd^g.) We assume that > nm to ensure that fiN possesses a continuous n-point 
density. Since is Gaussian, the density D^i^{b, a; z) is Gaussian as well, and it is described 
by the covariance matrix 



(57) A, 
where 



AN tdN 



nk \ jdN* /^N 
nk nk 



(5J 



A^,= [^s,{z^)s,,{zv')) 
Bl={'E^s,{z^)§^] 



(59) 



C^n^.= (E||(.-)|5(.^') 

jj' = l,...,k; = l,...,n; g, g' = 1, . . . , m. 

By (IID and (|5|), 

^nfc = {Sjj'SN{z^,zP')) , Sn{z,w) = i 1 + ^ZrW; 

\ r=l 

Bl = [5,,,SnA^p, y )) , SnA^^ = Nz,. ( 1 + 5^ 2;, 
Cl = {5,rSN,A^'. z^')) , Sm,A^. w) = N{N - \)W,z,, I 1 + ^ 

/ m 

+5,,. AT 1 + ^ 



^ N 



r=l 

rWr I , 



r=l 

N-1 



r=l 



The n-point zero correlation functions A'^ for the SU(m+ l)-polynomial /c-tuples Sn can 
be computed by substituting (^) into formulas (^OD and (|42|) . (Alternately, we can compute 
the zero correlation functions with respect to the Euclidean volume on C™" by setting 7^ =Id 
in (E2|).) 



Remark: Note that the one-point correlation function, or the zero-density function, is 
constant, since it is invariant with respect to the group SU(m + 1). Indeed, by Bezout's 



24 PAVEL BLEHER, BERNARD SHIFFMAN, AND STEVE ZELDITCH 

theorem and 



(60) |Z,|(l)=Vol(K)= / 7^a;"^-'= = fi2™-2fcdegK = ^^ 

JVs 



2m-2k 



where 

(61) n2i = Vo\cr = ^. 



Hence, 

(62) K^^iz 



We can also use our formulas to compute directly: By (p9D, 

(63) < = [s,Ai + \^\Y) , B^, = [s,,NzAi + \z\Y-') 



[6,,,N[{N - i)z,z,> + (1 + + kr)"^"') • 

Hence by (|g), 

(64) Af, = (5,yiV[(l + - -zM{^ + l^n^-^) = (5,yiV(l + \z\'rg,A^)) . 

In the hypersurface case (fc = 1), we compute 



= I / V 

11 7r(l + IzPW \ 



aiqlqq'O'lq 

q,q' = l 



N ^ Nm 
77- ^ 



TT — ' TT 



as expected. For > 1, we have A^(0) = NI where / is the unit matrix, and (^2]) yields 

^1.(0) = ^ ^det (^^a,,a,,, j ^ ^ j = _ , 
which agrees with (|62D. 



3. Universality and scaling 

Our goal is to derive scaling limits of the n-point correlations between the zeros of random 
fc-tuples of sections of powers of a positive line bundle over a complex manifold. We expect 
the scaling limits to exist and to be universal in the sense that they should depend only on the 
dimensions of the algebraic variety of zeros and the manifold. Our plan is the following. We 
first describe scaling in the Heisenberg model, which we use to provide the universal scaling 
limit for the Szego kernel (Theorem p.l|) . Together with Theorem |2.4| , this demonstrates the 
universality of the scaling-limit zero correlation in the case of powers of any positive line 
bundle on any complex manifold. 
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3.1. Scaling of the Szego kernel in the Heisenberg group. Our model for scaling is 
the Szego kernel for the reduced Heisenberg group described in §1.3.2| . Recall that for the 
simply-connected Heisenberg group H™', the scaling operators (or Heisenberg dilations) 

are automorphisms of H™ ( [[b'o|] p^]). Since the Szego kernel H of H"* is the unique self- 
adjoint holomorphic reproducing kernel, it follows that it must be invariant (up to a multiple) 
under these automorphisms. In fact, one has (|0, p. 538]): 

(65) U{6rX,6ry) = r-'"'-'U{x,y) 

The condition for a dilation 6r to descend to the quotient group H^^ is that r^Z C Z, or 
equivalently, r = a/ZV with N G Z+ . Note however that 5 ^ is not an automorphism of H™^ 

and there is no well-defined dilation by . 

The scaling identity ( p3| ) descends to in the form 

(66) U^{x,y) = N^Uf{5^x,6^y) 
with 

(67) Hf (x, y) = i^e^(^-^)e^^(-'^)e-il^-l^ , x = iz,9),y = {w, ^) . 

(Recall (|2T|).) Informally, we may say that the scaling limit of H^ equals H^. Since scaling 

by is not well-defined on H™^ it is more correct to say that H^ is the scaling of 

the scaling limit kernel. 

3.2. Scaling limit of a general Szego kernel. We now show that H^ is the scaling limit 
of the A^-th Szego kernel H^v of an arbitrary positive line bundle L — > M in the sense of the 
following "near-diagonal asymptotic estimate for the Szego kernel." 

Theorem 3.1. Let zq & M and choose local coordinates in a neighborhood of zq so that 
0/i(^o) = ^ d.Zj. Then 



iV-™H;v(2o + ^,^;^o + ^,^) = ^e*('^-^^+*^^"-"^-tl"-"l^ +0(Ar-i/2) 



uY{u,e;v,ip) + 0{N- 



-1/2N 



To prove Theorem |3.1| , we need to recall the Boutet de Monvel-Sjostrand parametrix 
construction: 

Theorem 3.2. PS| , Th. 1.5 and §2.c] Let H(x, y) be the Szego kernel of the boundary X of 
a bounded strictly pseudo-convex domain Q in a complex manifold L. Then: there exists a 
symbol s G ^"(X x X x IR+) of the type 

oo 

s{x,y,t) ~ ^t'^-''sk{x,y) 

SO that 

poo 

U{x,y) = / e''^^''^y\s{x,y,t)dt 
Jo 

where the phase ip G C°°(X x X) is determined by the following properties: 
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• x) = ^p{x) where p is the defining function of X. 

• dxip and dyip vanish to infinite order along the diagonal. 

The integral is defined as a complex oscillatory integral and is regularized by taking the 
principal value (see [ [BS|| ). The phase is determined only up to a function which vanishes to 
infinite order aX x = y and its Taylor expansion at the diagonal is given by 



(68) ip{x + M, x + f ) = - 



[x] 



The Szego kernels IItv are Fourier coefficients of 11 and hence may be expressed as: 
(69) UN{x,y)= / e-'^^e'''^^''''^yh{rex,y,t)dedt 



JO JO 

where vg denotes the action on X. Changing variables t \—>- Nt gives 

/»oo /»27r 

(70) UN{x,y) = N / e'^^--'^+'^^''<>''^y^^s{rgx,y,tN)dedt. 

Jo Jo 



We now fix zq and consider the asymptotics of 

= N I I e^^(-^+*^(-+^'^^^°+7F'0)),(;,o + ^,^;^o + ^,0),tiV)rf^^rft 



Un{zo + 0; Zq + 0) 



JO 



In our setting the phase takes the following concrete form: We let h{z, w) be the almost 
analytic function on M x M satisfying h(z,z) = WeiUh'^iz)- The function h{z,w) is defined 
by 

, ^d''+^h{z,z) . .u'^v^ 

We consider the complex manifold Y = L* and we let {z, A) denote the coordinates of 
^ e F given by ^ = A(e^)^. In the associated coordinates {x,y) = {z, X,w,fi) on F x y, we 
have: 

(73) p{z, X) = 1 — h{z, z)\X\'^ , tplz, X,w, p) = -{1 — h{z,w)Xfi) . 

i 

We consider Q = {p < 0} and X = dQ = {p = 0} . We may assume without loss of generality 
that h{z,w) = h{w,z) since h{z,z) is real so we could replace h by ^h{z,w) + ^h{w,z). On 

X we have h{z,z)\X\'^ = 1 so we may write A = h(z,z)~^e'^'^, and similarly for p. So for 
(x, y) = {z, (p, w , ip') & X X X we have 



(74) ilj(z,(f,w,f') = ^ 



1 _ h{z,w) 

^Jh{z,z)^Jh{w,w) 
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It follows that 



U V 



(75) 



1 



J9 



We now assume that is a normal frame centered at zq. By definition, this means that 
(76) h{zo) = 1, dhl, = dhl, = 0. 

We furthermore assume that our coordinates {zj} are chosen so that the Levi form of h is 
the identity at zq: 



(77) 



"(^0; Zq) — 5a(3 ■ 



dz'^dz'^ 

(This is equivalent to specifying that uj{zq) = | J2j dzj A dzj.) Then by (|7 



(7; 



h{zQ + 



Now let us return to the phase. It is given by 



id 



(79) t 1 ^ ^-^ -e 

hi^o + ^r^o + ^)-^Hzo + ^rzo + ^)-^ 

By ([78D, the phase (^) has the form: 



-ie. 



(80) 



U ■ V \u\ 

2' ' 



2' ' 



^3/2^ 



It is now evident that Hn^Zq + -^,0;zo + :^,0) is given by an oscillatory integral with 

phase {t[l — e*^] — i6); the latter two terms can be absorbed into the amplitude. 

Thus we have: 
(81) 

Un(zo + ^,0-zo + ^,0) 

= N j;^ J^"" e.^(t[i-e»«]-^e)giM-i|.p-i|.P]+o(7v-V2)^^^^ ^ ^. ^ 0), tiV)rf^(it. 

We may then evaluate the integral asymptotically by the stationary phase method as in ||Ze|| . 
The phase is precisely the same as occurs in Iliy{x,x), and as discussed in the single 
critical point occurs at t = 1, ^ = 0. We may also Taylor-expand the amplitude to determine 
its contribution to the asymptote. Precisely as in the calculation of the stationary phase 
expansion in |Ze|, we get: 



82 



^nIzq + 



u 



0;^o + 



,0) 



Finally, we note that 



1 



u ■ V 



\u\ \v\ 

II 2' I 



u-v-^\u\^-^\v\^ + 0(A^"~^). 



n — t>p + i'^iu ■ v) 
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which completes the proof of Theorem |3.1| . □ 

3.3. Universality of the scaling limit of correlations of zeros. We are now ready to 
pass to the scaling limit as ^ oo of the correlation functions of sections of powers 
of our line bundle. To explain this notion, let us consider the case k = m where the zeros 
are (almost surely) discrete. An m-tuple of sections of will have N"^ times as many 
zeros as m-tuples of sections of L. Hence we must expand our neighborhood (or contract 
our "yardstick") by a factor of N"^^'^. Let z° G M and choose a coordinate neighborhood 
U E M with coordinates {zj} for which z° = and u;(z°) = | Tlq^^q ^ ^^q- We define the 
n-point scaling limit zero correlation function 

K^kmiz) = Jim j^K^, , z={z\...,z-)e (C^n . 



We show below (Theorem ^.4| ) that this limit exists and that -ft'^m is universal by passing 
to the limit in Theorem ^.4| , using Theorem First, we need the following fact: 



Lemma 3.3. Let 2;^, . . . , be distinct points of C™. Then 

det (nf (z^ 0; z^', 0)) = e"^ l^^l' det (e^"-'"') ^ . 
Proof. We consider the first Szego projector on the reduced Heisenberg group 

(83) nf : £'(H™d) ^ ^i(H^d) ~ jC^C"', e-l^l') H 0(C"^) , 
where 

/:2(c-,e-l^l') = {/ G CLiCn ■■ j^M?e-\'\"dz < +00} . 

(See the remark at the end of § |1.3.2| .) Its kernel can be written in the form 

00 

(84) nf(z,^;ti;,V.) =e^(^-^) J]/„(^)^H, 

where the fa form a complete orthonormal basis for £^(C"^, e"'^'^) fl (9(C™). (E.g., {/„} 
can be taken to be the set of monomials [cj-^...j^zl^ ■ ■ ■ z-f^Y ia^ct, li^{z, 0; w, 0)) is just 
a "weighted Bergman kernel" on C™.) We now mimic the proof of Lemma p.3| , except this 
time we have an infinite sum over the index a; this sum converges uniformly on bounded sets 
in X C™ since the sup norm over a bounded set is dominated by the Gaussian-weighted 
norm (by the same argument as in the case of the ordinary Bergman kernel on a bounded 
domain). We then obtain a nonzero vector {vi, . . . ,Vn) G C™ such that J2p'^pfa{z^) = for 
all a. But then Ylp'^pfi^^) ~ polynomials / on C™, a contradiction. □ 

We can now show the universality of the scaling limit of the zero correlation functions: 

Theorem 3.4. Let {L, h) be a positive Hermitian line bundle on an m- dimensional compact 
complex manifold M with Kdhler form uo = |6/ij let S = H^{M,L^)^ (k > I), and give S 
the standard Gaussian measure fi. Then 



J^nk^nk I y]Y'---' J ~ ^nkm{^ ' \^/N 
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where K!^i^^{z^, . . . , z^) is given by a universal rational function in the quantities z^, z^, e^^'^^ , 
and the error term has i^'^ order derivatives < on each compact subset S C (C'^)„, for 
all£> 0. 

Proof. By taking the scaling limit of (0), we obtain 

^ ^ p„,„(nf (^^ y ), <nf (^^ z^'), <nf (^^ y ), « n?(^^ z^')) 



^5) 



k: 



nkm \ 



TT 



kn 



det (nH(,P,^P')) 



fe(n+l) 



l<p,p'<n 

Indeed, since the coefficients of A„ are either of degree 1 in the coefficients of C„ or of degree 

and 



2 in the coefficients of -B„, we see by the proof of Theorem 2.4, using 



Theorem |3.1| , that the leading term of the asymptotic expansion of K^f^ is N""^ times the 
right side of 



nk 

The bound on the error term follows from Theorem |3.1| and Lemma 



k: 



Substituting into (|85|) the values of n^(zP, z^') and its horizontal derivatives obtained from 
i^) (with N = 1) and ( |67D and canceling out common factors of e^l^''!^/^ and vr, we obtain 



nkm\ 



TX 



kn 



det (e^^-^"^') 



l<p,p'<n 



fc(n+l) 



(86) 



nkm 5 -^g 5 ^ ) 



^kn [det (e^-^-')]'=("+i)' 
where Qnkm is a universal polynomial (homogeneous of degree k{n+l) in each of the variables 



and with integer coefficients). 



□ 



Remark: As we remarked previously, formula (^) is valid for any connection, so we can 
replace the left invariant vector fields with their right-invariant counterparts to obtain 



^7) 



vr 



kn 



det (e^^-^^') 



q ' 1 

-1 fc(n+l) 



l<p,p'<n 



4. Formulas for the scaling limit zero correlation function 

We now apply the formulas from §§2.2- p^ to transform ( ^71) into explicit formulas for 
^nkm- We use the right-invariant connection so that df^ = Z^. Indeed, by the proofs of 



Theorems |2.4| and p.4| (which use formulas (^OD , (^21), (^6D-(^)), formula ( pTf ) becomes 

1 / " 



TT^"- det Ankm , -, 
\P=1 



(aPan 



where 

(89) ^nkm = Cnkm ~ ^nkm^nkm^nkm 
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with 

(m 
r=l 

Bnkm = {6jfSqf{zP,zP')) , Sq'{z,w) = Zgf exp ^ ZrW^ , 



,r=l ^ 

m 



(90) 



,r=l 



,r=l 



Cnfcm = {6jj'Sqq'{zP, Z^')) , (z, w) = {6qqf + WqZqi) CXp | ^ 

= l,...,k; p,p' = l,...,n; g, g' = 1, . . . , m. 

The metric tensor in (H3) becomes a unit tensor in the scahng hmit, so there is no 7*^ on 
the right in (B3). 

Because n{^is invariant with respect to unitary transformations and equivariant with 
respect to translations (i.e., n^(2 + u,w + u) = e*^^^""''e~*^''"'""^n^(z, w)), the scahng hmit 
zero correlation K!j^^ is invariant with respect to the group of isometric transformations — 
unitary transformations and translations — of C". 

In particular, the limit one-point zero correlation, or the zero-density function, is constant, 
since it is invariant under translation. Indeed by (0), ^ia,™ = e'^' Ik and Ai^m = e'""' Ikm, 
where J^, resp. Ikm, denotes the unit k x k, resp. (km) x [km), matrix. Thus by (^) and 
the Wick formula, 



(91) Kr,Jz) 



m! 



^k^k\z\^ yZ-.'-^^ I / 7r^(m-fc)!' 



Thus we define the normalized n-point scaling limit zero correlation function 
(92) K^^^{z) = [K'^^) ^K^p.^{z) = f — -j K^i^^{z) . 



Remark: These formulas also follow from § p.4| . For example, equation (^Tj) is a consequence 
of ( |62D since 

K^kmi^) = ^K^ki^) . 
Furthermore, using the notation of §2.4| , we observe that 

TV 



I = lim ( 1 + A^"^ V Zrw; 1 



lim Sjsf ( -A=, — ^ 1 = lim I 1 + ^ "N " ZrWr 1 = Siz, w) 



y^^^N-'S^qq, = Sqq,{zM . 



Equations ( ^3]) provide an alternate derivation of (|90D. 



UNIVERSALITY AND SCALING OF CORRELATIONS BETWEEN ZEROS 



31 



4.1. Decay of correlations. Explicit formulas for the correlation functions can be 



obtained from (pSf), (|90|) and the Wick formula. We shall illustrate these computations for 



the cases n = 2, A; = 1, 2 in 5S|4.2H4.3| below. We now note that the limit correlations are 



"short range" in the following sense: 

Theorem 4.1. The correlation functions satisfy the estimate 

K^f^.^{z^, . . . , z"-) = 1 + 0{r^e~''^) as r ^ oo , r = min |z 



Proof. We use formula (p5|), which comes from (|88D -(|89D as in the proof of Theorem 3^. To 
determine the matrices A, B, C, we let = Z^, = (instead of the right-invariant 
vector fields we used above). Recalling (p2|), we have: 



(94) = S,A4'-4K'^ 



By (§7p, 



1 p = p' 



-P' - I 0(e-^ /2) p^p' ' 
B = 0(re-"'/2)^ 
C = / + 0(rV^'/2)^ qp. 



Recalling (]40|), we have 



(95) A = / + Oir^e-'"/^) , AH = 1 + 0{r\-'') . 



We now apply formula (^Bp; note that the Wick formula involves terms that are products 
of diagonal elements of A, and products that contain at least two off-diagonal elements of 
A. The former terms are of the form 1 + 0(r^e~'' ), and the latter are 0(r^e~^ ). Similarly, 
det A = 1 + 0{r^e~^ ). The desired estimate then follows from (^2]). □ 



We shall see from our computations of the pair correlation below that Theorem O! is 
sharp. The theorem can be extended to estimates of the connected correlation functions 
(called also truncated correlation functions, cluster functions, or cumulants), as follows. The 



n-point connected correlation function is defined as (see, e.g., [|GJ| , p. 286]) 



I 

(96) f-^(^\ ...,zn = $^(-i)'+^(/ - 1)! n ^n%(^^ • • • , ^"0, 

G j=l 
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where the sum is taken over all partitions G = {Gi, . . . ,Gi) of the set (1, . . . , n) and Gj = 
(zi, . . . , in^). In particular, recalling that = 1, 

= ^3kmi^ 5 ^ 5 ~ ^2km{^ 5 ^ ) ~ ^2km{^ 5 "~ ^2km{^ ? + ^ , 

and so on. The inverse of (^) is 

z 

(97) n = Y.X{ TuM'^'' ' • • • ' ) 

G j=l 

(Moebius' theorem). The advantage of the connected correlation functions is that they go 
to zero if at least one of the distances |z* — goes to infinity (see Corollary |4.3| below). In 



our case the connected correlation functions can be estimated as follows. Define 
(98) d{z\ ...,zn= maxj] k^^'^ - ^/W I^e-I^'^^'-^'^'IV^. 

l&L 

where the maximum is taken over all oriented connected graphs Q = (V, L, i, /) such that 
V = {z^, . . . , z"') and for every vertex z^ G V there exist at least two edges emanating from 
z^ . Here V denotes the set of vertices of Q, L the set of edges, and and /(/) stand for 
the initial and final vertices of the edge /, respectively. Observe that the maximum in ( p8|) 
is achieved at some graph Q, because te~*/^ < 2/e < 1 and therefore the product in ( pSD is 
less or equal (2/e)l-^l which goes to zero as \L\ oo. 

Theorem 4.2. The connected correlation functions satisfy the estimate 



T„^„(z\ . . . , z") = 0{d{z\ . . . , z^) as max \z^ - 



oo 



P,Q 



provided that miUp^g |z^ — > c > 0. 

This theorem implies Theorem |4.1| because of the inversion formula (|97D. To prove the 
theorem let us remark that we can rewrite (88) (using the Wick theorem) as a sum over 
Feynman diagrams. Namely, for the normalized correlation functions K!^^{z^, . . . ,2;") we 
have that 

(99) K-jz\ ...,.") = [^"1"///"^'^'^ y: ^^(-^ 

Get ^nkm ^ 

where the sum is taken over all graphs JF = (V, L, i, /) (Feynman diagrams) such that 
V = (2;^, . . . , z^) and the edges / G L connect the paired variables a*^^'*, oj/g/* in a given term 

of the Wick sum for K^j^{z^, . . . , z^). The function Ajr[z^^ . . . , 2;") is a sum over all terms 
in the Wick sum with a fixed Feynman diagram J-'. In other words, to get Ajr[z^, . . . , z^) 
we fix pairings (a^^, oJ^ig/) prescribed by JF and sum up in the Wick formula over all indices 
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j, q at every a^. A remarkable property of the connected correlation functions is that they 
are represented by the sum over connected Feynman diagrams (see, e.g., [pJ|| ), 



rpOO / 

nkm\ 



detA 



nkm 



Since det Ankm > Ci > and |Ajy_j,| < C2 < +00 when min 



z'^l > c > 0, we conclude 



from (|40|), (Q) and (0) that for all connected Feynman diagrams JF, 
(101) A^(^\...,^") = 0(rf), 

where d = (i(2;^, . . . , z"") is defined in (0). Summing up over JF, we prove Theorem O. □ 
Corollary 4.3. The connected correlation functions satisfy the estimate 

as i? ^ 00 , R 



z")=0(i?V«/^) 
z1 > c> 0. 



max \z^ — z'^l 
p,<i 



provided that miup^g \zP 

4.2. Hypersurface pair correlation. We now give an explicit formula [( |110| )] for pair 
correlations in codimension 1 {k = 1, n = 2). The case m = 1 of this formula coincides, 
as it must, with the formula given by ||Han|| and ||BBL|| for the universal scaling limit pair 
correlation for SU(2) polynomials. In another paper |[BSZ|| , we gave a different proof of ( |110| ) 
using the Poincare-Lelong formula. 

Since the scaling-limit pair correlation function K^^lz^jz"^) is invariant with respect to 
the group of isometrics of C™', it depends only on the distance r = — so we can set 
z^ = and z"^ = (r, 0, . . . , 0). To simplify notation, we shall henceforth write A = A2km., B = 

B2km, C = C2km, A = A2km- 



In this case, ( pOf ) reduces to 
(102) 



'-ill' 

B = (5,\); 

The matrix 
(103) 

is given by 
(104) Al\ - 





r re^ 

1 1 

1 (l + r^)e^ 



(BP, 








[CI?,) =5, 



qq' 



1 1 

1 



g,g'>2. 



A = [Al\,) = C- B*A-'B 



/ e" - 1 - u e" - 1 - ue" \ 

e" - 1 e« - 1 

e" — 1 — we" e^" — e" — ue" 



\ e" - 1 e" - 1 / 

where u = r'^ = — z'^\'^ . By (|88|), (p2D and the formula for A in ( |102| ), we have 



(105) 



-^21m(^ ? ^ ) 
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By the Wick formula (see for example, |S], (1.13)]) 
(106) 



m2(e" — 1) 
1 



.9=1 



g,q'=l 



.9=1 



W' = l 



m?{e'^ — 1) 

Substituting the values of A^f^, given by (|104|) , we obtain 
(107) 



1 



l-U 



m?{e^ — 1) 

'e" — 1 — ue" 



+ 



e" - 1 

2 



+ m — 1 



g,g'=l 



g2?i _ g« _ 



+ (m - l)e" 



- 1 



+ (m - 1) 



I 1 2|2 
U = \z — Z \ 



M2(e2« + e") - 2M(e2" - e") + m2(e" - 1)^6" + m(e" - 1)^ 



After simplification, 

(108) K^,m{z\z\ 
Putting u = 2t and writing 

(109) K^,^{z\z') = kU\^'-z'\), 
we then obtain 

, , [^(m^ + m) sinh^t + t^l cosht — (m + l)tsinht (m — 1) 

(110) Ki„(r) = ^ ^—^ ^ + ^ ^, „ 

sinh t 2m 2 

The case m = 1 of formula (|110|) was obtained by Bogomolny-Bohigas-Leboeuf [ [BBL|] and 



t 



Hannay |Han 



As r — i> oo, 



(111) 



— 2(m2 + l)r2 + m{3m + 1) 



e-' + 0{r^e 



The following expansion of the correlation function was obtained from (|110|) using Maple 



TM. 



m — 1 _^ m — 1 1 (m + 2) (m + 1) 1 (m + 4) (m + 3) 



2m 



-t 



2m 6 



t 

90 



m^ 



_^ 1 (m + 6)(m + 5)^5 1 (m + 8)(m + 7)^7 



945 
1 



9450 



93555 



(m + 10) (m + 9) ^9 _ 

2 



691 (m+ 12) (m + 11)^11 



638512875 



+ 



2 (m + 14) (m+ 13)^13 



18243225 m2 



In particular, in the one-dimensional case we have 

2 36 720 16800 



(112) 



1 



^14 + 
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4.3. Pair correlation in higher codimension. Next we compute the two-point correla- 
tion functions for the case k = 2. For A; > 1, we have 

(113) A = {A%,) = {6,yA;,), fi = (5iW) = fe'5p%'), c = {qr^%,) = {s,,,c;:'^,) , 



where A^,, 5^,^,, C^y are given by (|102D . It follows that 
(114) A = (Af;,,,) = {6,,A;f^,) 

where A^^ is given by ( |104| ). 



By (HP, 
(115) 



7r2^(e" - 1)^= 



det 



a]a], 



det 



j,j'=l,...,k 



2 9 



j,j'=l,...,k 



V 



9=1 



where m = = — as before. Observe that the random variables a^^ and a^,^, are 
independent if either j ^ j' or g 7^ g'. 
RecaUing (p^), we write 



(116) K^,^{z\z') = kU\z'-z'\). 

When k = 2, (|115D reduces to the following 

{alal){alal) - {a\al){alal) (afaDialal) - {afaDialaf) 



(117) 



By the Wick formula, 
(118) 



m'^{'m — l)2(e" — 1)^ 

dii — (^21 — di2 + 

im?{m — l)2(e" — 1) 



2 ' 



where 



(119) 



similarly. 



d 



11 



'a\a\){a\al){alal){alal) 



5^ 



E Al^Al^A^^Ai +2 5:A-A»^^ 



(120) 



(121) 



d 



12 



>ia})(a2«2)(«i«2)(«2«i) 



^21 



[Ag]^) (e, A}^) + 2 (E, AgA?^A^^) (E.AlVj + [A^^A^ 

(aia2)(«2«i)(«i«i)(a2a2 

[A^^) (E, Ag)' + 2 (E, A^A^^) (e, Ag) + E, [A^^A^ 
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,i7T^ /'^l7T^ /'^2~2\ /^s"! ^ 



(122) 



d22= {ia^al){a2al){a^al){a2aj) 



I A 2 

[A^ ') (E, [Ag] V 2 A}^AgA?^A^^ + (e, A^A^^) . 



Substituting the values of the matrix elements of A we then obtain 



[m 



2 



m)e^'' + 2(m - l)e" + 2 Aue''[{m - l)e" + l](m + 1) 



(123) 



(e" — l)^m(m — 1) (e" — l)^(m — l)m^ 



2M2e«[(m - l)e2" + 2me" + 1] 



\ ; TTT ^ — n ■, u = r 



2 



(e" - l)4(m - l)m 



2 



As r ^ oo, 



, , 2[r^-2(m + l)r2 + m(m + l)le-^ 
124 n2m = l + ^ ^ ^ + O r^e 



As r 0, 
(125) 



m-2 4 m-2 , Sm^ - 7m + 12 (m - 2)(m + 2)(m + 1) , 

K2m{r) = r-^ + ^ u + - 7^7 TtH 

m m 12[m — l)m 12(m — Ijm^ 

(m + 3)(m + 2) 4 (m - 2)(m + 4)(m + 3) g 
^ 240(m-l)m ^ 720(m - l)m2 



When m = 2 the asymptotics reduce to 
(126) ^22(r) = ^ + ^- 



4 24 288 4800 96768 
and in this case K22 is a series in r^. 
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